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Abstract 

When an inextensible elastic rod is ‘injected’ through a sliding sleeve against a fixed constraint, 
conhgurational forces are developed, deeply influencing the mechanical response. This effect, 
which is a consequence of the change in length of the portion of the rod included between the 
sliding sleeve and the hxed constraint, is theoretically demonstrated (via integration of the 
elastica) and experimentally validated on a proof-of-concept structure (displaying an interesting 
force reversal in the load/deflection diagram), to provide conclusive evidence to mechanical 
phenomena relevant in several technologies, including guide wire for artery catheterization, or 
wellbore insertion of a steel pipe. 
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1 Introduction 

The buckling of a piece of paper hitting an obstacle when ejected from a printer, the insertion of a 
catheter in to an arteryjB] or of steel piping in to a wellbore |151 1201 F6] . even the so-called ‘inverse 

spaghetti problem ’ miniiB] are all examples of mechanical settings where an elastic rod is forced 
out of a constraint and pushed against an obstacle. During this process the length of the rod subject 
to deformation changes, and when the rod is injectec0 through a sliding sleeve, a conhgurational or 
Eshelby-like force is generated, that until now has been ignored in the above-mentioned problems 
[IHllIT]. This force has recently been evidenced for a cantilever beam by Bigoni et al. |3]. It has 
also been shown to influence stability |2], used to design a new kind of elastically deformable scale 
[5] and to produce a form of torsional locomotion jl]. 

‘Corresponding author: Davide Bigoni, bigoni@unitn.it; +39 0461 282507 

^ Similarly to a fluid in a vein, ‘injection’ means that a portion of an elastic rod is forced to enter in a fixed space 
between two constraints. 

^The concept of configurational force was developed by Eshelby [^I10| in the context of solid mechanics to motivate 
the possibility of motion of defects, such as dislocations, cracks and voids. In structural mechanics such a force was 
demonstrated for the first time in [3]. 
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The aim of this article is to provide direct theoretical and experimental evidence that the effect 
of conhgurational forces on the injection of an elastic rod is dominant and cannot be neglected. It 
leads, in the structure that will be analyzed, to a force reversal that otherwise would not exist. The 
considered setup is an inextensible elastic rod, clamped at one end and injected from the other 
through a sliding sleeve (via an axial load), so that the rod has to buckle to deflect (Fig. left). 



Figure 1: Left: Injection of an elastic planar rod through a sliding sleeve via an axial load P. The rod of total 
length I has bending stiffness B, is clamped at its lower end and inserted into a frictionless sliding sleeve 
(at a distance L from the clamp) at the upper end. Due to symmetry, the system can be analyzed as 
four equal cantilever beams of length ieq/4 loaded with a dead load P + /2B, where /2B is the 

axial configurational force generated at the sliding sleeve. Right: the experimental set up realized with an 
elastic rod (1), load cells (2), sliding sleeve (3) and movable crosshead (4). 


The elastic system is analytically solved in Section]^ through integration of the elastica mw- 
It is shown that ‘Eshelby-like’ forces strongly influence the loading path and yield a surprising force 
reversal, so that certain equilibrium configurations are possible if and only if the applied force changes 
its sign. The change of sign is shown to occur when the rotation at the inflexion points exceeds 7r/2 
(corresponding to Al/L ~ 1-19), a purely geometric condition independent of the bending stiffness 
B and distance L. Furthermore, it is also proven that during loading two points of the rod come 
into contact and, starting from this situation (again defined by a purely geometric condition), the 
subsequent configurations are all self-intersecting elastica. All theoretical findings have been found 
to tightly match the experimental results presented in Section and obtained on a structural model 
(designed and realized at the Instability Lab of the University of Trento http://ssmg.unitn.it), see 
also the supporting electronic material. 
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2 The elastica and theoretical behaviour of the structure 


The structure shown in Fig. (left) is considered axially loaded by the force P. The curvilinear 
coordinate s G [O, is introduced, where I is the total length of the rod, together with the rotation 
held 9{s) of the rod’s axis. Denoting with a prime the derivative with respect to s, the total potential 
energy of the system can be written as 


^out 


v{eis),iont) = 


B 


O'{s) 


1 2 


^As-P 


I — / cos0(s)ds 


( 1 ) 


where /out £ [L, 1) is the length of the deformed elastic planar rod comprised between the two 
constraints (the clamp and the sliding sleeve) and B is the bending stiffness of the rod. 

With reference to e <C 1 (Bigoni et al. |2) the rotation and the length of the rod comprised 
between the two constraints can be written as the sum of the helds evaluated at equilibrium, 0eq 


and /eq, and their variation 0var and 


^(•^T) — ^eq(®) T C^var('S), /out(f^) — /eq T C/\i 


( 2 ) 


so that consideration of the boundary conditions at the sliding sleeve, 0eq(/eq) = 0 and 0(/out) = 0, 
leads to the compatibility equation 


^var(/eq) — ^eq(^eq) ^varj 

from which the hrst variation of the functional V follows as 


(3) 


6,V = - 


// 1 B t 

Be'' + Psin6'eq(s) 6»var(s)ds - —O'^Meqfl 


(4) 


0 


The fact that the distance L between the two constraints remains constant for every deformation is 
expressed by the constraint 

^out 


L= / cos0(s)ds, 


(5) 


SO that an additional compatibility condition is obtained 

^eq 

/var = / Sin6'eq(s)6'var(s)ds, 


( 6 ) 


and the vanishing of the hrst variation S^V (Q leads to the governing equation of the elastica 
^eq(s) + (p + ■^^eq(^eq)^'j sin0eq(s) =0, S G [0, /eq]- 


(7) 
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The bending moment M is proportional to the rod’s curvature, M{s) = so that the equilibrium 

condition Q reveals the presence of an axial conhgurational force proportional to the square of the 
bending moment: M(Zeq)^/2-B, see (Bigoni et al. j^). 

Restricting attention to the hrst buckling mode and exploiting symmetry, the investigation of 
equation Q can be limited to the initial quarter part of the rod, s G [0, ^eq/4], which may be treated 
as the cantilever rod shown in the inselt of Fig. The equilibrium conhguration can be expressed 
as the rotation held 0eq('S), solution of the following non-linear second-order differential problem 


Cq(s) +P^sin6'eq(s) = 0, 
0eq(O) = 0, 



s G 




( 8 ) 


where the parameter p, representing the dimensionless axial thrust, has been introduced as 




P . oUo? 


Integration of the elastica and the change of variable 


V = sm-. 


• \ ■ ^eq('S) 

V sm (p[sj = sm —^— 


(9) 


( 10 ) 


where 0eq (^eq/4) = 0, yields the following differential problem for the auxiliary held (p 


4> (s) = p\J 1 — sin^ (p{s), 

</( 0 ) = 0 , 



s G 




( 11 ) 


Further integration of 01 leads to the relation between the load parameter p and the rotation 
measured at the free end of the cantilever 6 (which is an inflection point )[^ 


p leq = 4% (u), (13) 

where % (v) is the complete elliptic integral of the hrst kind. 

The configurational force, included in p, is a function of the curvature at the ends of the rod 
comprised between the two constraints, namely 0eq(O) = ^eq(^eq)) which can be obtained, through 

^Would the configurational force /2B be neglected, the following solution is obtained 



( 12 ) 
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a multiplication of equation (j8j)i by 6^^ and its integration, together with the boundary condition 
^eq (^eq/4) = ^, aS 

<q(0) = ^J2p^{l-cose) = 2pv, (14) 

so that equation may be rewritten as 


I = ohl - 2‘>") . 


The rotation held 0eq{s) can be expressed through the inversion of the relation (10)2 as 

0eq(s) = 2arcsin (usn {sp,v )), 


(15) 


(16) 


while the axial and transverse equations describing the shape of the elastica are obtained from 
integration of the following displacement helds 


xi{s) = / cos0eq(r)dr, 

Jo 


X 2 {s) = / sin6»eq(r)dr, 


as 


xi(s) = —s H— { E [am {sp, v), u]} 
P 

X 2 {s) = — [1 -cn{sp)] , 

P 


(17) 


(18) 


where the functions am, cn, and sn denote respectively the Jacobi amplitnde, Jacobi cosine amplitude 
and Jacobi sine amplitude functions, while E{x,v) is the incomplete elliptic integral of the second 
kind of modulus v. Note that equations (16) and (18) are valid for the entire structure, s G [0; ^eq] 
Although the problem under consideration seems to be fully determined by equations (0) 
and (18), the length /eq is unknown because it changes after rod’s buckling. This difficulty can be 
bypassed taking advantage of symmetry, because the axial coordinate of the rod, for every unknown 
cantilever’s length /eq/4, is xi(/eq/4) = T/4, so that equation (18)i gives 




am 


eq 


P,V ,V 


L 

4’ 


(19) 


and therefore, inserting equation (13) in equation (19), it is possible to obtain the relation between 
the load parameter p and the angle of rotation at the free edge of the cantilever 9 as 


p = y {2 E [am (X(v),v ), u] — lK(r;)} 

1j 


( 20 ) 


The applied thrust P, normalized through division by the Eulerian critical load of the structure 
P„ = 47r2p/L2^ is calculated from equation ( |l5[ ) as a function of the kinematic parameter 9 through 
equation (|10[)i as 

P 4 (1-2^2) 




TT^ 


{2 E [am (JC(i;), u), u] — J<C(u)} . 


( 21 ) 


Fig. shows load P (divided by Per) versus the rotation at the inflection point 9. Also shown, for 
comparison, the dashed line representing the structural response when the conhgurational force is 
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neglected, equation (12|. It is noted that, although the critical load is not affected by the presence 
of the configurational force (because the bending moment is null before buckling), the behaviour is 
strongly affected by it, so that the unstable postcritical path exhibits a force reversal when 9 > 90° 
(confirmed also by experiments, see Sect. |^, absent when the configurational force is neglected. 

According to equation (21|, a self-equilibrated, but deformed, configuration (P = 0) exists for 
9 = 90°, where the rod is loaded only through the configurational force M'^/2B. 



90 117.54 130.71 

0 20 40 60 80 100 120 140 




Figure 2: Loading path of the structure sketched in Fig. (left), expressed as the dimensionless applied dead load 
P/PcT versus the rotation 9 at the inflexion point, equation (151, shown with a red line. The loading path 
is compared with that obtained by neglecting the configurational force (dashed black line), equation (12l. 
Two deformed equilibrium configurations have been sketched in the inset showing that a compressive 
thrust (A) or tensile load (B) are necessary to guarantee equilibrium for 9 < 90° and 9 > 90°, respectively. 
For both solutions, after the condition 9 « 117.54° (marked with a green star) is met, a self-intersecting 
elastica occurs, while load reversal does not occur when the configurational force is neglected. The rotation 
value 9 ~ 130.71° is a limit condition corresponding to an infinite amount of injected rod length, A1 —>■ oo, 
and a null applied load, P = 0. 


The length A/, measuring the amount of rod ‘injected’ through the sliding sleeve 

Al = (eq - L, (22) 


can be easily computed from equations (13) and (20) in the following dimensionless form 

Al 


JCiv) 


L 2 E [am (lK(t;), u), u] — lK(u) 


- 1 , 


(23) 


whereas, considering the relationship (14), the analytical expression for the ‘Eshelby-like’ force 
(divided by Per) becomes 


__ = ^ {2 E [am , u] - %{v)y 


(24) 
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When the rotation at the inflexion points exceeds 9 ~ 117.54° (corresponding to A/ ~ 5.61L), the 
elastica self-intersects, so that the presented solution holds true for rods capable of self-intersecting 

(as shown in mm)- 

Finally, it is noted that the condition 


2 E [am (iK(u), u), u] — %{v) = 0 


(25) 


occurs for 9 ~ 130.71°, representing a limit condition for which the injected rod length becomes 
inhnite, Al —?• oo, and the applied load becomes null, P = 0, as evident from equations (231 and 
(211, respectively. 


3 Experimental 

The structural system shown in Fig. (right) was designed and manufactured at the Instabilities 
Lab (http://ssmg.unitn.it/) of the University of Trento, in such a way as to be loaded at prescribed 
displacement Al, with a continuous measure of the force P. 

The displacement was imposed on the system with a loading machine MIDI 10 (from Messphysik). 
During the test the applied axial force P was measured using a MT1041 load cell (R.C. 500N) 
and the displacement Al by using the displacement transducer mounted on the testing machine. 
Data were acquired with a NI compactRIO system interfaced with Labview 2013 (from National 
Instruments). The elastic rods employed during experiments were realized in solid polycarbonate 
strips (white 2099 Makrolon UV from Bayer, elastic modulus 2350 MPa), with dimensions 650 mm 
X 24 mm x 2.9 mm. The sliding sleeve, 285 mm in length, was realized with 14 pairs of rollers from 
Misumi Europe (Press-Eit Straight Type, 20 mm in diameter and 25 mm in length), modihed to 
reduce friction. The sliding sleeve was hxed to the two columns of the load frame, in the way shown 
in Fig. (right). 

The influence of the curvature of the rollers on the amount of the conhgurational force generated 
at the end of the sliding sleeve was rigorously quantihed in |3] and found to be negligible when 
compared to the perfect sliding sleeve model. 

A sequence of six photos taken during the experiment with L = 360 mm is reported in Fig. 
together with the theoretical elastica, which are found to be nicely superimposed to the experimental 
deformations. 

Experimental results are reported in Fig. and compared with the presented theoretical solution 
and the solution obtained by neglecting the conhgurational force (dashed line). While the former 
solution is in excellent agreement with experimental results performed for two different distances 
between the two constraints (L = 360 mm and L = 410 mm), the latter solution reveals that 
neglecting the conhgurational force introduces a serious error. Note that two different lengths L 
have been tested only to verify the robustness of the model, because size effects are theoretically 
not found with the assumed parametrization, which makes results independent also on the bending 
stiffness B. 

The dimensionless applied load P/Pcr versus dimensionless length Al/L (Fig. left), measuring 
the amount of elastic rod ‘injected’ into the sliding sleeve, conhrms a load reversal at Al/L ~ 1.19. 
We note that the conhgurational force does not inhuence the qualitative shape of the elastica, but the 
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Figure 3: A sequence of photos of experiments performed on the structure sketched in Fig. (right) in the undeformed 
(A) and deformed configurations (B-F). The six photos show the displacement-controlled experiment at AZ = 
{0; 0.14; 0.41; 0.82; 1.21; 1.58} L, corresponding to measured loads P = {0; 0.62; 0.26; 0.06; 0;—0.02} Per. 
The calculated postbuckling deformed configurations are superimposed to the photos as dashed lines. 


amount of deflection, so that neglecting it leads to a strong and completely unacceptable estimation 
of the amplitudes. 

Finally, the dimensionless conhgurational force M^/2i?Pcr reported as a function of Al/L (Fig. 
on the right) shows its strong influence on postcritical behaviour, such that this force grows to make 
up one third of the rod’s critical load. 

4 Conclusion 

Instabilities occurring during injection of an elastic rod through a sliding sleeve were investigated, 
showing the presence of a strong conhgurational force. This force causes a force reversal during a 
softening post-critical response and has been theoretically determined and experimentally validated, so 
that it is now ready for exploitation in the design of compliant mechanisms. Moreover, conhgurational 
forces as investigated in this paper can also be generated through growth or swelling, so that the 
can play an imortant role in the description of deformation processes of soft matter. 
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